AD  A 03 1 577 


DNA  3967Z 


DATA  REDUCTION  AND  ANALYSIS- 
EXPERIMENTAL  IEMP/SGEMP  ELECTRON 
BEAM  STUDY,  FIRST  INTERM  REPORT: 

40  KEV  BEAM  ENERGY 

Kaman  Sciences  Corporation 
P.O.  Box  7463 


Colorado  Springs,  Colorado  80933 
December  1975 
Interim  Report 

CONTRACT  No.  DNA  001-75-C-0167 

APPROVED  POR  PUBLIC  RELEASE; 
DISTRIBUTION  UNLIMITED. 


THIS  WORK  SPONSORED  BY  THE  DEFENSE  NUCLEAR  AGENCY 
UNDER  RDT&E  RMSS  CODE  B323075464  R99QAXEB08953  H2590D. 


Prepared  for 
Director 

DEFENSE  NUCLEAR  AGENCY 
Washington,  D.  C.  20305 


UNCLASSIFIED 


security  CLASSIFICATION  OF  this  page  (When  Dele  Entered) 


REPORT  DOCUMENTATION  PAGE 


2.  GOVT  ACCESSION  NO. 


4 TITLE  (end  SuhNlI#) 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3.  RECIPIENT'S  CATALOG  NUMBER 


5.  TYPE  OF  REPORT  & PERIOD  CQVEREO 


J? ATA  REDUCTION  AND  £N  AL  Y SIS -Jj^tP  ER I MENTAL 
IE^g/SGEMP  ELECTRON  £EAM  JTUDY,  JIRST 
jNTER~:mwfaSP0RT : 4^  KEV  BEAM  JJNERGY^ 


AnT  NUMBERf S) 


PERFORMING  ORGANIZATION  name  ANO  adoress 

Kaman  Sciences  Corporation  / 

P.0.  Box  7463  y 

Colorado  Springs,  Colorado  80933 


11.  CONTROLLING  office  NAME  ANO  ADORESS 

Director 

Defense  Nuclear  Agency 
Washington,  D.C.  203Q5 


SS (It  different  from  Controlling  Of lice) 


is.  NUMBER  OF  PAGES 
60 


15.  SECURITY  Class  (ot  thla  report ) 

unclassified 


15j.  DECLASSI  FI  CATION/ DOWNGRADING 
schedule 


16-  DISTRIBUTION  STATEMENT  (Ol  thla  Report) 


Approved  for  public  release;  distribution  unlimited. 


IT  DISTRIBUTION  STATEMENT  (oi  the  abmtract  entered  In  Block  20,  If  different  from  Report) 


1#  SUPPLEMENTARY  NOTES 

This  work  sponsored  by  the  Defense  Nuclear  Agency  under  RDT&E 
RMSS  Code  B323075464  R99QAXEB08953  H2590D, 


IS.  KEY  WORDS  (Continue  on  reveree  aide  If  necaaamry  and  Identify  by  block  number) 

SGEMP  Simulation 
Electron  Beam  Experiments 
IEMP  Simulation 
Electromagnetic  Pulse  (EMP) 


20  IMi^TRACT  (Continue  on  r#v*ri*  aide  If  nec eatery  and  identify  by  block  number) 

J*This  reporters  Kaman  Sciences'  first  formal  interim  report 

Y>which  presents  reduced  data  from  th^  series  of  internal  electro- 
magnetic  pulse  (IEMP)  simulation  experiments^fcarried  out  by  the 
Ballistic  Missile  Defense  System  Command,  eady  and  it 

also  presents  reduced  time  history  and  tabular  data  on  selected 
pressures  for  the  nominal  40  keV  mean-energy  injected  beam. 


. 

••  '^r 
■%% 


- V-'  ‘ fx'- 

■":m0' 

_ ,-^r;V£^v ; 


;v  ,^v'  * 
••-#  • 


UNCLASSIFIED 

SECUBITV  CLASSIC  1CATIQN  Of  THIS  t»AGE(>Wl«n  Q«l«  gnUfttf) 

20.  ABSTRACT  (Continued) 

Details  of  the  experimental  apparatus  and  techniques  are 
found  in  the  following  reference:  "A  Study  of  IDMP/SGEMP 

Phenomenology  Using  Electron  Beam  Simulation;"  also,  additional 
discussion  is  presented  in  Appendix  B in  this  report.  Addi- 
tional similar  interim  reports  will  cover  the  remaining  two 
sets  of  data  comprising  the  nominal  60  keV  mean-energy  beam 
into  the  18  cm  length  chamber  and  nominal  125  keV  beam  into  the 
30  cm  length  chamber. 

j 

A final  report  will  contain  analyses  and  correlation  of  the] 
t cee  data  sets  covered  in  the  interim  reports  as  well  as  the 
60  keV  data  set  previously  reduced  and  reported  on  in  the 
above-mentioned  reference.  The  final  report  will  also  contain 
results  of  computer  calculations  using  Kaman  Sciences'  two- 
dimensional  code  EBBTIDE  for  comparison  with  appropriate 
determined  quantities. 


accession  nr 

NT  It 

White  Seetlrt 

DEC 

toll  Sectfis 

□ 

UNANK5UWE0 

□ 

JUSTIFICATION 

oiSTntaaTm/AVAiutsiuTY  «oo 

o»*T  Ttaii.  wi/yj tfltiiT 


PREFACE 


The  bulk  of  the  data  processing  techniques  and  codes 
employed  in  reducing  the  data  presented  in  this  report  were 
developed  under  the  leadership  of  Dr.  J.  S.  Duval  prior  to 
his  leaving  Raman  Sciences  to  join  the  U.  S.  Geological  Survey. 
Additional  discussions  with  him  during  the  course  of  the  work 
presented  herein  are  also  acknowledged  with  pleasure. 


1 


TABLE  OF  CONTENTS 


e 


! 


Page 


PREFACE 

I.  INTRODUCTION 7 

II.  CHARACTERIZATION  OF  THE  INJECTED 

ELECTRON  BEAM 9 

III.  PRESENTATION  OF  DATA 13 

Magnetic  Fields  13 

Transmitted  Current . 14 

Collected  Wall  Currents  22 

Radial  Current  Densities 26 

IV.  PEAK  VALUES 26 

V.  ERROR  ANALYSIS 32 

General  Considerations 32 

Conversion  of  Oscilloscope  Traces 

to  Volts  vs  Time 32 

Transmitted  Current  Measurements 

(Total  and  Radially  Resolved)  ....  34 

Wall  Segment  Currents  35 

Magnetic  Fields  37 

Flush-Mounted  Current  Probes 40 

APPENDIX  A 4 3 

APPENDIX  B 47 


3 

l — — 


x. 

"7  • 


• r v/.va*  % 




;, k}-  • ‘ - - 


• • r- ^4-.  •;.?•' 

v fm ' 


- 

- .0# 


■ w-,. 

■* 


Figure 

1 


2 

3 


[LLUS THAT IONS 
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upper  and  lower  curves  are  standard 
deviations 

"Typical"  machine  current  monitor  output. 
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in  Figure  1 
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DATA  REDUCTION  AND  ANALYSIS- 
EXPERIMENTAL  IEMP/SGEMP  ELECTRON  BEAM  STUDY, 
FIRST  INTERIM  REPORT:  40  keV  BEAM  ENERGY 


I.  INTRODUCTION 

The  following  report  is  Kaman  Sciences*  first  formal 
interim  report  to  the  Defense  Nuclear  Agency  on  work  sup- 
ported under  contract  number  DNA001-75-C-0167 . Presented 
herein  are  reduced  data  from  the  series  of  internal  electro- 
magnetic pulse  (jEMP)  simulation  experiments  carried  out 
under  sponsorship  of  the  Ballistic  Missile  Defense  System 
Command,  contract  number  DASG60-74-C-0015 , early  in  1974. 

The  experiments  carried  out  were  designed  to  explore 
under  controlled  laboratory  conditions  the  phenomenology 
governing  the  behavior  of  a uniformly  distributed  beam  of 
electrons  injected  axially  into  a test  chamber  of  cylindrical 
geometry  and  to  provide  a reliable  body  of  experimental  data 
against  which  to  evaluate  existing  or  evolving  system 
generated  electromagnetic  pulse  (SGEMP)  computational  codes. 

The  cylindrical  experimental  chamber  was  instrumented 
in  such  a manner  as  to  provide  correlated  time  history  data 
on  transmitted  currents  as  a function  of  radial  location  and 
wall  currents,  radial  current  densities  and  magnetic  fields 
as  a function  of  axial  location.  Additionally,  radial 
electric  field  measurements  were  attempted  using  short 
monopole  stubs,  but  plasma  sheathing  complications  rendered 
those  data  unusable. 
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A 30  cm  diameter  monoenergetic  electron  beam  generated 
using  the  Simulation  Physics  Incorporated  SPI  5000  machine 
was  injected  axially  into  the  experimental  chamber  through  an 
aluminized  mylar  membrane.  Suitably  gasketed  to  complete  the 
vacuum  tight  integrity  of  the  chamber,  this  transmitting  ancde 
of  the  machine  diode  formed  the  entrance  plate  and  allowed  the 
drift  chamber  to  be  mainuained  at  a fixed  pressure  for  each 
shot.  The  pressure  range  explored  extended  from  less  than 
one  micron  to,  in  some  instances,  twenty  torr.  The  time 
histories  of  the  quantities  enumerated  above  along  with 
electron  beam  machine  diagnostic  voltage  and  current  outputs 
were  monitored  for  each  shot  on  twenty-three  separate  data 
channel."  (each  equalized  to  350  MHz  bandwidth)  and  displayed 
on  Tektronix  519  oscilloscopes.  The  total  body  of  data 
acquired  corresponded  to  three  different  mean  energies  of  the 
injected  beam  -40,  60,  and  125  keV;  two  different  chamber 
lengths,  IB  and  30  cm;  and  the  pressure  range  stated  above. 


Details  of  the  experimental  apparatus  and  techniques  are 
to  be  found  in  reference  (1).  Additional  discussion  is  also 
presented  in  Appendix  B of  this  report.  This  report  will 
present  reduced  time  history  and  tabular  data  on  selected 
pressures  for  the  nominal  40  keV  mean-energy  injected  beam. 

The  length  of  the  test  chamber  was  30  cm  and  the  injected 
current  pulse  length  was  30  nsec  FWHM.  Additional  similar 
interim  reports  will  cover  the  remaining  two  sets  of  data 
comprising  the  nominal  60  keV  mean-energy  beam  into  the  18  cm 
length  chamber  and  nominal  125  keV  beam  into  the  30  cm  length 
chamber.  A final  report  will  contain  analyses  and  correlations 
of  the  three  data  sets  covered  in  the  interim  reports  as  well 
as  the  60  keV  data  set  previously  reduced  and  reported  on  in 
reference  (1),  The  final  report  will  also  contain  results  of 
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computer  calculations  using  Kaman  Sciences'  two-dimensional  code 
EBBTIDE  for  comparison  wtih  appropriate  experimentally  determined 
quantities.  EBBTIDE  provides  a self -consistent , fully  electro- 
magnetic solution  to  charged  particle  eleourodynamic  problems. 
Comparisons  of  EBBTIDE  calculations  to  date  have  been  against 
electron  beam  experimental  data  corresponding  to  end-on  photon 
illumination  of  a right  circular  cylinder  and  against  underground 
test  data  of  a right  circular  cylinder  enclosing  a sphere  located 
on  the  axis  of  the  cylinder.  Results  from  the  three-dimensional 
code  TRIFIC,  which  uses  the  same  basic  algorithms  as  EBBTIDE, 
have  also  been  compared  against  underground  test  data  corre- 
sponding to  side-on  illumination  of  concentric  right  circular 
cylinders.  In  all  cases,  agreement  has  been  to  within  a factor 
of  two  in  amplitude  and  waveshape. 

The  remainder  of  this  report  has  been  organized  into  four 
main  sections  contained  in  Sections  II,  III,  IV  and  v.  Section 
II  describes  the  process  which  has  been  used  to  obtain  the 
characterization  of  the  injected  electron  beam  mean  energy 
and  current  (as  a function  of  time) . Section  IT  contains 
time  history  waveforms  of  various  data  categories.  Section 
IV  contains  summary  tables  and  plots  of  selected  data  from 
Section  III.  Finally,  Section  V discusses  the  experimental  and 
processing  errors  inherent  in  the  data  of  Section  III. 

I 1 . CHARACTERIZATION  OF  THE  INJECTED  ELECTRON  BEAM 

A reliable  working  understanding  of  the  nature  of  the 
injected  beam  is  of  course  crucial  in  interpreting  and  applying 
the  experimental  results  to  actual  systems  situations. 
Additionally,  an  accurate  specification  of  the  injected  beam 
in  terms  of  energy,  current  density  and  angular  distribution 
is  required  input  to  any  computer  code  which  is  intended  to 
simulate  IEMP  phenomena  occurring  within  the  chamber.  A 
detailed  comp. .ison  of  experimental  vs  computed  results  using 
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Haitian's  EBBTIDE  code  is  to  be  presented  in  the  final  report  on 
this  project.  The  curves  presented  in  this  section  will  serve 
to  characterize  the  injected  beam  in  sufficient  detail  to  allow 
useful  interpretation  of  the  experimental  results  presented 
within  the  report. 

As  is  discussed  in  reference  (1)*,  the  injected  60  keV 
electron  beam  used  in  generating  the  data  presented  therein 
can  be  considered  quite  uniform  over  its  entire  diameter. 

Recent  evidence,  however,  has  appeared  indicating  that  under 
the  operating  conditions  of  the  SPI  PULSE  5000  electron  beam 
machine  used  to  generate  shorter  pulse,  lower  energy  data  such 
as  that  presented  herein,  the  injected  beam  may  assume  a degree 
of  nonuniformity.  This  nonuniformity  is  manifested  in  a current 
density  distribution  in  which  the  beam  center  is  less  dense 
than  the  beam  edge  by  something  on  the  order  of  a factor  of 
one  half.  This  circumstance  is  of  critical  importance  in  the 
context  of  code  comparison  against  the  experimental  data  and 
will  have  to  be  considered  when  defining  the  code  input. 

Three  additional  quantities  beyond  uniformity  must  be 
specified  in  order  to  fully  characterize  the  injected  beam 
which  drives  the  experiment;  beam  energy,  current  density  and 
angular  distribution  of  the  injected  electrons.  Appendix  III 
of  reference  (1)  discusses  in  detail  the  considerations  and 
procedures  involved  in  converting  the  machine  diagnostic  out- 
puts corresponding  tc  shank  voltage  and  total  diode  current 
to  the  more  pertinent  quantities  of  interest,  namely,  injected 
beam  energy  and  current  density.  Briefly,  what,  jr  required  is 
to  correct  the  machine  voltage  monitor  output  f jc  the  inductive 
voltage  drop  between  the  monitor  itself  and  the  diode  gap  in 


* Reference  (1).  Duval,  J.S.,  Rich,  W.F. , and  Clarke,  C.A.  "A 

Study  of  IEMP/SGEMP  Phenomenology  Using  Electron  Beam  Simulation.’1 
Kaman  Sciences  Report  (K-74-126U (R) , November,  1974. 
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which  the  electron  beam  is  generated*  With  the  correction 
applied,  the  resultant  voltage  waveform  represents  the  time 
history  of  the  accelerating  potential  applied  to  the  electrons 
emitted  from  the  diode  cathode.  At  each  instant  in  time,  the 
electrons  impinging  upon  the  anode  plane  have  an  essentially 
monoenergetic  distribution  which  corresponds  to  the  potential 
difference  that  has  accelerated  them  across  the  diode  gap. 

In  order  to  arrive  finally  at  the  energy  distribution  of  the 
injected  beam,  the  impinging  electron  beam  must  be  transported 
through  the  .25  mil  aluminized  mylar  diaphram  and  wire  mesh 
which  constitutes  the  diode  gap  anode.  Reference  (1)  contains 
curves  generated  by  the  Monte  Carlo  code  EPIC  which,  when 
applied  to  the  indicant  electron  beam,  yield  mean  energy, 
angular  distribution  and  current  density  of  the  beam  which  is 
ultimately  injected  into  the  test  chamber. 

Figure  1 presents  the  "typical"  machine  voltage  monitor 
waveform  associated  with  the  group  of  data  presented  in  this 
report.  This  waveform  was  constructed  by  computing  the  mean 
and  standard  deviation  values  of  the  individual  machine  voltage 
waveforms  (calibration  factor  of  9.8  kv/v)  at  two  nanosecond 
intervals  for  each  shot  in  the  series  included  in  this  report. 
The  center  curve  represents  the  mean  or  the  "typical"  wave- 
form, while  the  upper  and  lower  curves  represent  the  ± a values. 

The  inductive  voltage  correction  described  above  depends 
upon  an  L -jJ-  term  in  which  I is  the  current  flowing  in  the 
diode  circuit,  namely  the  current  measured  on  the  machine 
current  monitor  (calibration  factor  of  64  amp/volt) . A 
"typical"  machine  current  waveform  was  generated  in  the  same 
manner  as  described  above  for  machine  voltage.  Figure  2 pre- 
sents the  resulting  curves.  The  mean  curve  shown  in  this 
figure  was  implemented  in  correcting  the  machine  voltage  in 
Figure  1.  The  curve  in  Figure  3 represents  the  resulting 
diode-gap  or  electron-accelerating  voltage. 
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Figure  4 presents  mean  injected- beam  energy  vs  time 
corresponding  to  the  "typical*'  waveform  in  Figure  1.  This 
figure  can  be  used  to  determine  the  mean  energy  of  the 
injected  beam  corresponding  to  the  machine  voltage  output  in 
Figure  1 for  any  instant  in  time. 

The  second  beam  quantity  of  interest  is  current  density. 
Figure  5 displays  the  injected  beam  current  measured  on  a 
solid  aluminum  plate  1 cm  behind  the  machine  anode  plane. 

(A  discussion  of  the  errors  in  the  injected  current  measure- 
ment is  to  be  found  on  page  172  of  reference  (1) . Recently 
available  data  on  the  injected  current  suggests  that,  due  to 
space  charge  limiting  very  close  to  the  injection  plane. 

Figure  5 may  be  in  error  on  the  low  side.) 

III.  PRESENTATION  OF  DATA 

Data  acquired  using  the  40  keV  beam  described  in  the 
previous  section  included  time- correlated  and  dependent 
moebius-lo op  B data,  total  and  spatially- resolved  transmitted 
currents,  collected  wall  currents  and  localized  radial  current 
density.  In  this  section  are  presented  time  history  curves  of 
data  from  each  of  these  categories.  It  should  be  noted  that 
during  processing,  the  polarities  of  some  of  the  sets  of  time 
history  data  presented  in  this  section  have  been  inverted. 

This  is,  of  course,  incidental  and,  once  noted,  should  not 
cause  difficulty  in  interpretation  of  the  data. 

Magnetic  Fields 

Because  of  the  relatively  low  signal-to-noise  ratio 
inherent  in  the  40  keV  data  considerable  difficulty  was 
encountered  in  the  processing  of  the  B loop  outputs  in 
this  data  group.  In  particular,  the  integration  process 
required  in  converting  the  moebius  loop  output  to  magnetic 
fields  is  quite  sensitive  to  systematic  digitizing  errors 
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resulting  from  imprecision  in  establishing  the  zero  reference 
as  discussed  in  Section  V.  This  is  of  course  true  for  all 
cases,  but  the  problem  is  compounded  from  a relative  error 
point  of  view  when  the  signal  amplitudes  are  small  as  in  this 
instance.  This  difficulty  is  discussed  in  Section  V. 

Figures  6a  through  6j  display  magnetic  fields  obtained 
over  the  pressure  range  of  5 to  800  microns.  In  these  and 
subsequent  figures,  the  individual  curves  are  displaced 
vertically  simply  as  a convenience  for  presentation  purposes 
so  that  the  progression  of  behavior  as  a function  of  pressure 
is  apparent. 

Transmitted  Current 

The  curves  in  Figures  7a  and  7b  display  the  total  trans- 
mitted currents  collected  on  the  solid  chamber  endplate  as  a 
function  of  time  and  cavity  pressure. 

Figures  8a  through  8m  display  the  currents  detected  on 
each  of  the  individual  concentric  rings  which  constitute  the 
chamber  endplate.  These  rings  permit  resolution  of  the  radial 
distribution  of  the  transmitted  current  into  five  regions. 
However,  the  rings  can  be  driven  in  an  antenna  mode  at  times 
in  the  presence  of  high  frequency  radial  electric  fields. 

(See  reference  (1),  p.  183.)  Hence  in  some  cases  they  do  not 
respond  as  simple  current  collectors  and  thus  some  care  is 
indicated  in  the  interpretation  of  their  output.  Table  I 
gives  the  dimensions  of  the  endplate  collection  rings,  which 
are  approximately  equal  in  area. 
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TWmSMITTED  CURRENT  (KA1 


RING  NO. 


IlHESSflkBSinUMS  s 

INNER 

OUTER 

- 

6.7 

6.8 

9.5 

9.6 

11.7 

11.8 

13.6 

13.7 

15.0 

Collected  Wall  Currents 

The  segmented-chamber-body  construction  of  the  test  chamber 
described  in  reference  (1)  allows  the  monitoring  of  total  current 
flowing  in  the  chamber  wall  at  five  successive  locations  down  the 
length  of  the  cavity.  The  subtraction  of  the  signal  developed 
across  the  resistive  gasket  material  in  one  section  from  that 
immediately  nearer  the  injection  end  of  the  chamber  yields  the 
net  current  collected  on  the  wall  segment  in  question.  Thus  the 
average  radial  currents  over  6 cm  segments  are  sampled  at  five 
locations  down  the  length  of  the  chamber.  Table  II  summarizes 
the  locations  of  the  various  sections.  These  results  are 
presented  in  Figures  9a  through  9 j , 


TABLE  II.  Locations  of  Wall  Segments.  Z = 0 
represents  the  injection  plane. 


Wall  Channels 

Segment  No.  Subtracted 


22A  - 15A 
15A  - 16A 
16A  - 17A 
17A  - 18A 
18A  - 19A 


Iz  Location  (cm)  ! 

From 

To 

0 

6 

6 

12 

12 

18 

18 

24 

24 

30 

22 


10.0 


Figure  9).  tfall  current*  collected  on  segmented' 

ri9ure  9i.  Well  current*  collected  on  **g»*nted—  chamber  alemnts. 

cheeber  elements. 


Radial  Current  Densities 


An  additional  set  of  measurements  were  performed  in  the 
40  keV  sequence.  These  consisted  of  sampling  the  radial  cur- 
rents at  seven  localized  points  by  means  of  flush-mounted 
current  probes  (FMCP)  mounted  in  the  chamber  wall.  Since 
these  probes  were  employed  only  in  the  40  keV  shot  sequence, 
a discussion  of  the  probe  construction  and  response  is 
included  in  this  report  as  an  Appendix.  The  response  analysis 
discussed  therein  indicates  that  the  FMCP  is  an  excellent 
device  for  sampling  the  impinging  current  densities.  Figures 
10a  through  lOj  present  data  obtained  using  these  sensors. 

We  have  elected  to  display  these  data  on  the  same  time  scales 
as  the  wall  segment  currents  to  facilitate  comparison  of  their 
wave  shapes. 

IV.  PEAK  VALUES 

Figure  11  displays  peak  magnetic  field  values  measured 
at  Z = 3,  15  and  27  cm,  that  is,  near  the  entrance,  middle 
and  end  of  the  drift  chamber.  These  data  are  plotted  as  a 
function  of  pressure. 

Figure  12  similarly  presents  the  peak  values  of  total 
transmitted  current  against  pressure.  In  Figure  13  are 
curves  of  collected  wall  current  peak  values  for  chamber 
sections  near  the  front,  middle  and  rear  of  the  test 
chamber . 


RADIAL  Ct'IUUMT  DdUIfV  (A/C* 


■uoiu  cumewt  »nom  <Vc**> 


Figure  10i.  Radial  currant  density  measured  by  the 
f luah-a*>unted  currant  probes  (PMC P) . 


8 

s 


rigure  10).  Radial  current  danelty  measured  by  the 
f lush«mounted  currant  probes  (PWCP) . 


figure  11.  Peek  magnetic  fields  a*  pressure 
near  front,  middle  and  rear  of 
chamber- 


^■2SS*3Uii x 


Table  III  summarizes  the  transmitted  current  distri- 
butions as  obtained  from  the  peak  values  of  collected  current 
on  each  of  the  individual  endplate  rings. 


TABLE  III.  Peak  currents  on  endplate 
collector  rings. 


Pressure 

(u) 

1 

Peak  Current 

(kA) 

f 1 ' 

Ring  1 

Ring  2 

Ring  3 

Ring  4 

Ring  5 

5 

.03 

.02 

.04 

.04 

.03 

10 

.05 

.04 

.05 

.05 

.05 

15 

.03 

.04 

.04 

.06 

.06 

20 

.07 

.06 

.07 

.07 

.07 

25 

* 

.06 

.05 

.06 

.05 

30 

.12 

.09 

.10 

.09 

.09 

35 

.15 

.13 

.13 

.09 

.12 

45 

.56 

.39 

.36 

.30 

.28 

55 

.52 

.41 

.44 

.35 

.28 

65 

1.  31 

1.16 

1.19 

.88 

.79 

200 

2.54 

1.92 

1.67 

.91 

.73 

400 

2.31 

1.60 

1.39 

.98 

.74 

800 

1.81 

1.28 

1.12 

.76 

.58 

* Bad  data  film. 
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V. 


ERROR  ANALYSIS 


General  Considerations 

Sources  of  error  have  been  divided  into  three  categories: 
experimental  error,  trace  digitization  error  and  propagation  of 
these  two  error  types  in  subsequent  processing.  The  third  area 
of  concern  deals  with  the  sensitivity  of  the  final  results  to 
inaccuracies  in  the  input  data. 

Conversion  of  Oscilloscope  Traces  to  Volts  vs  Time 

In  order  to  set  the  time  scale  on  each  scope,  a sine 
wave  of  known  frequency  was  recorded  on  each  channel  photograph 
prior  to  each  shot.  The  distance  (from  an  arbitrary  zero  point) 
to  each  successive  peak  of  this  waveform  was  digitized.  From 
these  data  the  mean  y'  and  standard  deviation  a"  of  the  distances 
xi+1~x^,  i=l,  2,  ...,n  between  successive  peaks  were  computed. 

Let  X be  the  random  variable  associated  with  distance 
errors  in  the  digitized  x (time)  axis.  We  assume  that  X has 
mean  0 and  variance  a2.  Let  y be  the  "actual"  sine  wave 
inter-peak  distance,  then  X-^  will  be  distributed  as  X + iy. 

Hence  Xi+x  - Xi  will  have  mean  u and  variance  2a2  given 
independence  of  X^  and  X^+i  . 

* 

The  statistic  y"  was  used  to  estimate  y and 
was  used  to  estimate  a.  c'  was  found  to  be  less 
than  .036  mm  on  nearly  every  shot.  Other  sources  of  error, 
such  as  non-uniformity  of  the  sine  wave,  have  been  ignored. 

These  sources  would  cause  a slight  under-estimate  of  data 
quality  rather  than  an  over-estimate.  It  is  assumed  that 
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voltage  errors  (vertical  distances)  are  distributed  the  same 
as  time  errors  (horizontal  distances) . 


Voltage  values  are  calculated  from  the  raw  digitized 
data  as: 


V . (fldo) 

m 


a s 


where : 


d is  the  vertical  coordinate  of  the  trace  on  the 
film  (mm) 

dQ  is  the  vertical  coordinate  cooresponding  to  0 (nun) 

m is  an  optical  magnification  factor  associated  with 

the  digitizing  process 

a is  signal  attenuation,  and 

. . . / v 

s is  scope  sensitivity  I 


Relativising  voltage  errors  AV  to  the  peak  voltage  encountered 

I v I max ' one  obtains 


o^A(d-dQ)^  was  approximated  by  .036  mm  (see  above)  and  ld-d0lmax 
was  tabulated  for  a sample  of  curves  from  each  type  of  data. 
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In  all  curves  investigated 


a 


< 5.5% 


Other  values  assigned  were: 


Hence  o 


(, 


AV 


» 6%  which  represents  the  error  inherent  in 


•max, 


the  digitized  data  upon  which  all  further  processing  was  performed. 


Transmitted  Current  Measurements  (Total  and  Radially  Resolved) 

These  measurements  result  from  measuring  voltages  developed 
across  resistive  shunts  and  therefore  involve  application  of  Ohm’s 
Law  in  the  form  I = V/R.  Percentage  errors  associated  with  peak 
values  were  evaluated  according  to  the  following  procedures. 

Letting  AI,  AV  and  AR  denote  random  variables  correspond- 
ing to  current,  voltage  and  resistance  errors  respectively  with 
variances  a2  (AI) , o2( AV)  and  a2(AR),  one  obtains 


As  determined  in  the  preceeding  section,  ivi  ) 

\i  I max/ 

The  worst  case  for  the  channels  under  consideration 

— ] < 8%  for 
\ 1 max/ 


in  this  section  was  about  .05.  Hence 


the  total  transmitted  current  measurements.  Recent  considera- 
tions relative  to  interpretation  of  the  current  collection 
ring  data,  however,  has  indicated  that  the  response  of  these 
rings  is  not  well  understood.  The  antenna  mode  response  of 
the  concentric  ring  configuration  appears  to  be  greater  than 
originally  determined.  It  is  not  possible  at  this  time, 
therefore,  to  complete  the  error  analysis  of  those  sensors. 


Wall  Segment  Currents 

The  net  wall  current  in  each  of  the  five  sections  com- 
prising the  body  of  the  drift  chamber  is  computed  as 
V2(t)/r2  - vj(t)/ri  where  vj  (t)  and  vz (t)  are  the  digitized 

voltages  associated  with  two  successive  sections  and  r,  and  r2 

are  the  corresponding  gasket  resistances.  The  error  analysis 
here  is  complicated  by  the  fact  that  relative  timing  discrep- 
ancies between  the  two  voltage  waveforms  must  be  accounted  for 
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Timing  errors  in  a given  trace  are  assumed  to  be  of  two  kinds. 
The  first  is  a systematic  error.  That  is,  due  to  inaccurate 
determination  of  the  timing  zero  point,  the  entire  waveform  is 
shifted  in  time  with  respect  to  true  zero  and  hence  with 
respect  to  the  waveform  with  which  it  is  to  be  differenced. 

The  second  error  source  is  a random  inaccuracy  in  the 
determination  of  the  time  of  each  individual  data  point. 

It  was  judged  that  these  random  errors  were  quite  small 
leaving  the  systematic  error  of  the  first  type  dominant. 

To  analyse  the  effect  of  the  systematic  timing  errors, 

let  i i ( t ) and  i2  (t)  be  correct  current  values  whose  difference 

i2(t)  - ii (t)  is  desired.  Due  to  errors  of  various  sorts, 

i2'(t)  is  obtained  instead  of  i2 (t)  and  ij"(t)  instead  of 

ii(t).  Allow  a relative  time  shift  error  of  a = 2.5  nsecs  and 

assume  that  there  exist  time  shifts  Si  and  S2  such  that 

|SjJ<a  and  | ik  (t) -i'  (t+S^)  | <ek  for  k - 1,2  and  all  t of 

interest.  In  other  words,  assume  that  i^(t)  time-shifted  by  a 

small  amount  Sk  is  a good  approximation  to  ik(t).  Define 

*>.  (t)  = max 
|s|^a 

| S | < a.  Hence  bk(t)  is  the  maximum  amount  that  i'ft)  can 
vary  as  it  is  time-shifted  by  amounts  varying  from  -a  to  a. 

The  difference  between  the  calculated  value  and  the 
desired  result  i2(t)  - ii(t)  is  of  interest.  The  triangle 
inequality  yields: 

|ia (t)-ij (t)-(ia(t)-ii (t))  | 


| i^ (S+t) -i' (t) | then  | i£ (S+t) -i' (t) | <bk (t)  for 


i|ia(t)-ia(S2+t) | + | i£ (S2+t) -i a (t) | + | i l (t) -i T (S i+t) | + | if (S i+t) -i 1 (t) | 


<.  b2  (t)  + e2  + bi  (t)  + ei 
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The  bound  b^(t)  measures  the  local  variation  in 
i'(t)  and  will  vary  with  t.  In  simplest  terms,  small  time 
shifts  will  cause  dramatic  changes  in  i^(t)  where  (t)  is 
large  and  i^(t)  has  steep  slopes. 

In  a number  of  cases  including  high  and  low  pressures, 
b^(t)  was  computed  for  a =2.5  ns,  k = 1,  2 and  t corresponding 
to  the  maximum  of  |i£  (t)— (t)  j . bj(t)  +b2(t)  was  found  to  be 
less  than  Q.7  kA  in  nearly  every  case.  Assigning  standard 
deviations  of  0.4  kA  to  £i  and  £2  one  obtains  a standard 
deviation  of  0.9  kA  for  the  total  error.  The  estimate  is 
believed  to  be  conservative  in  the  sense  that  it  bounds  the 
true  standard  deviation  but  that  0.6  kA  is  probably  a more 
realistic  estimate. 

Magnetic  Fields 

The  magnetic  field  results  were  calculated  from  the 
moebius  loop  data  using  the  formula: 


B ( t ) 


L 

“ Tkr 


t 

v(t)+S  j v(x)dx 
o 


where: 

A = (1.08  ± .01)  x 10-“m2 

R = 25Q 

L = (16.9  ± .01)  x 10“3hy 

v(t)  is  the  digitized  voltage  trace  corrected  for  data 

channel  attenuations.  That  is,  v(t)  represents  the 
output  voltage  of  the  sensor  itself. 
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The  presence  of  errors  in  v(t)  can  be  handled  by 
replacing  v(t)  with  v(t)+e(t)  where  e(t)  is  an  error  term. 
Define 


B*  (t ) 


L 

Iar 


t 


v(t)+e  (t)+^ 


Hence,  |B*(t)-B(t)| 


L 

2AR 


t 


£(t)+S 


(T)dx 


_ L|e(t)  I 

* 2ar 


+ 


_l 

2A 


t 

J e(i)dT 

o 


L[e(t)  1 
2AR 


+ 


t 

2A 


max 

o<f<t 


|e{x)  | 


Because  errors  of  a small  level  shift  and/or  rotation  are 
judged  present  in  the  process  of  reading  the  photographic 
images,  one  must  allow  for  cumulative  error  rather  than 
assuming  that  the  error  will  "average  out"  over  time. 

This  cumulative  error  corresponds  to  the  term: 

t 

7K  j e(T}dT 

o 
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IS iiEZ: 


The  following  example  illustrates  the  problem  of 
cumulative  error.  Given  a uniform  level  shift  error  in  v(t) 
of  2%  of  its  peak  value,  a simple  computation  shows  that 


t 


_1_ 

2A 


J e(x)dr 


o 

could  accumulate  an  error  equal  to  the  peak  value  of 


Lv  (t ) 
2AR 


in  time 


34  ns.  At  t 


34  ns  it  would  not  be 


uncommon  to  find  values  of  B*(t)  of  roughly  this  same 
magnitude . 


In  summary,  the  values  calculated  for  B(t)  decrease  in 
accuracy  with  advancing  time  and  can  be  rather  unreliable  even 
at  times  as  early  as  35  ns.  The  values  for  B(t)  at  times  near 
zero  approach  the  accuracy  of  the  digitized  voltages  since  the 
errors  in  A,  R and  L are  small. 


What  this  means  in  practice  for  this  set  of  data  in 
which  both  the  signal  levels  are  small  and  pulse  durations  are 
short  is  that  in  those  magnetic  fields  for  which  the  peak  values 
occur  beyond  40  nsec,  the  peak  values  are  not  very  reliable. 

It  is  estimated  that  a bound  on  the  error  on  the  computed 
value  of  the  B fields  is  probably  between  50  and  75%  at  40 
nsec.  There  is  evidence,  however,  that  the  error  is  not 
nearly  that  bad  in  much  of  the  data;  but  quantifying  it  in 
each  case  is  extremely  difficult.  One  reason  for  believing  the 
error  is  much  less  is  that  on  many  shots  in  which  there  are 
three  loops  azimuthally  distributed  at  given  Z coordinate, 
the  resulting  calculated  B fields  agree  to  within  15%.  In 
other  cases  the  disagreement  is  marked,  and  a late  time  ramp 
attributable  to  a fixed  zero  offset  in  the  raw  data  is  apparent. 
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Attempts  to  simply  subtract  out  the  associated  zero  were  not 
satisfactory  and  it  was  felt  that  further  tinkering  with  the 
data  to  produce  better  "agreement"  was  scientifically 
unwarranted. 

It  should  be  noted  that  the  60  and  125  keV  data  sets 
have  much  less  error  associated  with  them.  This  is  due  to 
two  factors.  First  those  data  have  significantly  better 
signal-to-noise  characteristics  in  that  the  signals  are  of 
both  greater  amplitude  and  longer  pulse  duration.  Secondly, 
the  digitizing  of  those  data  was  performed  using  much  more 
accurate  equipment  and  processes  than  the  present  set. 


Flush-Mounted  Current  Probes 

The  formula  J (t)  = Kiv(t)  + K2~X.;—  was  used  to  cal- 
r dt 

culate  the  radial  current  densities  from  the  flush-mounted 
current  probe  voltage  data  with  K3  = 129.9  and  K2  = 158.4. 

(See  Appendix.)  The  derivative  at  t2  was  approximated  by 

evaluating  the  slope  of  a parabola  passing  through  (tj,v(ti)), 
(t2,v(t2))  and  (t 3 , v (t 3 ) ) where  ti,  t2  and  t3  are 
consecutive  time  points.  A new  parabola  was  fit  to  successive 
sets  of  three  points  advancing  along  the  data  points  one  at  a 
time . 


Here  the  sensitivity  of  the  derivative  calculation  to 
errors  in  the  input  data  is  of  over-riding  concern.  The 
result  Jr(t)  is  least  sensitive  to  errors  in  v(t)  when  either 

the  term  is  small  or  the  term  K23^(t)  is  at  least  small 
dt  dt 

with  respect  to  the  term  Kiv(t). 
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In  order  to  evaluate  the  sensitivity  of  the  derivative 
calculation  to  errors  in  the  input  times,  the  derivative  was 
written  as  a function  g (t 1 ,t2 , t 3 ,vi ,V2 ,v3)  and  the  gradient 
Vg  was  computed  for  actual  data  points.  For  the  most  part  the 
standard  deviations  of  the  timing  errors  were  approximately 
0.25  nsec.  It  was  found  that  times  near  the  steepest  part  of 
v(t)  could  yield  errors  exceeding  20%  of  the  resulting  0r(t) 
peak  value  even  though  the  time  error  (At 1 ,At2 ,At3)  satisfied 
the  requirement  ( (At  1 ) 2+ (At 2 ) * + (At 3 ) 2 ) % < 0.25  nsec.  The 
derivative  computation  is  a good  deal  less  sensitive  to  errors 
in  v(t)  where  v(t)  is  steep.  At  times  when  the  curve  v(t) 
flattens  out,  the  accuracy  of  the  Jr(t)  values  approaches  the 
accuracy  of  v(t). 

In  summary,  the  flush-mounted  current  probe  data  has  an 
error  bound  of  about  20%  but  is  probably  routinely  good  tc  10%. 

i 

c 


i 
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APPENDIX  A 


FLUSH-MOUNTED  CURRENT  PROBE 

The  flush-mounted  current  probes  (FMCP)  employed  in 
this  experiment  were  designed  to  measure  the  radial  electron 
current  density  impinging  on  localized  areas  of  the  test 
chamber  walls.  The  probes  themselves  consist  of  a current 
collection  element  which  was  a rod  of  diameter  7 mm  insul- 
ated from  its  surrounding  support  structure  by  a 1.05  mm 
layer  of  epoxy.  The  entire  probe  assembly  was  inserted 
into  the  wall  structure  of  the  chamber  so  that  the  face  of 
the  probe  was  flush  with  the  inside  wall  of  the  chamber  and 
the  support  structure  and  chamber  walls  were  common. 

Terminals  on  the  collection  element  and  the  support 
structure  provide  the  signal  path  to  a 50  ohm  terminated 
line  leading  directly  to  a Tektronix  519  oscilloscope. 

There  are  two  potential  modes  of  excitation  of  the 
sensor  which  have  been  analysed:  direct  current  collection 
of  the  radial  component  of  the  electron  current  impinging 
on  the  sensor  face  (the  desired  response)  and  the  antenna 
response  due  to  the  axial  component  of  the  electric  field 
in  the  vicinity  of  the  sensor.  An  additional  excitation  mode 
is  due  to  displacement  currents  in  the  vicinity  of  the  sensor 
but  this  mechanism  is  judged  to  be  negligible  and  was  not 
explicitly  considered.  An  analysis  of  the  probe  response  has 
been  performed  assuming  ramp  functions  of  appropriate  slopes 
for  the  E and  J driving  functions  based  on  the  equivalent 
circuit  model  for  the  probe  and  its  50  ohm  load  shown  in 
Figure  1-A. 
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Figure  1-A.  Equivalent  circuit  model  of  the 
flush  mounted  current  probe. 

In  this  Figure,  v(t)  represents  the  antenna  voltage, 

C the  antenna  capacitance,  i (t)  represents  the  collection 

a 

current  on  the  center  rod,  Cp  the  capacitance  formed  by  the 
current  collection  electrode  and  the  rest  of  the  support 
structure,  RL  is  the  50  ohm  load  and  VL  is  the  measured 
voltage  output  of  the  probe  system.  The  driving  ramp 
functions  were  chosen  to  have  slopes  similar  to  the  rise 
times  of  Jr  and  Ez  observed  in  the  experiment.  Table  1 
tabulates  the  analysis  results.  V_  is  the  response  voltage 

a 

at  the  load  due  to  the  antenna  effects  of  the  probe,  V is 
that  due  to  the  impinging  current  and  V is  the  total  response. 
The  table  indicates  that  for  times  as  early  as  1 nsec,  the 
error  in  attributing  the  entire  response  to  the  collected 
current  is  less  than  3%.  At  later  times  the  error  becomes 
totally  negligible.  Therefore,  the  equivalent  circuit  may 
be  simplified  even  further  to  that  shown  in  Figure  2-A  the 
circuit  equation  of  which  is 
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r APPENDIX  B 

INSTRUMENTATION  DETAILS 

In  this  section  some  additional  details  of  the  instru- 
mentation of  the  test  chamber  not  included  in  reference  1 are 
presented. 

Magnetic  Field  Sensors 

Standard  moebius  strip  loops  with  potted  gaps  were  employed 
as  the  magnetic  field  sensors.  A good  discussion  of  the  princi- 
ples of  operation  of  the  loops  themselves  can  be  found  in  the 
AFWL  EMP  Sensor  and  Simulation  Note  Number  7.  Figure  1-B  shows 
the  data  channel  configuration,  while  Figure  2-B  presents  the 
schematic  details  of  the  balun  employed.  Figure  3-B  is  a 
photograph  of  the  sensor  itself.  Figure  3 of  reference  1 shows 
the  loops  mounted  in  the  drift  chamber  walls. 

Current  Collection  Sensors 


Voltage  signals  proportional  to  current  through  the  resis- 
tive gasket  material  separating  the  individual  chamber  segments 
were  picked  off  directly  by  attaching  50Q  RG-58  cables  to  the 
chamber  segment  flanges.  Figure  4-B  illustrates  this.  The 
equal  area  current  rings  were  isolated  from  the  solid  endplate 
also  by  resistive  gaskets,  and  the  voltage  signals  developed 
by  collected  current  on  each  of  those  was  similarly  tapped. 
Resistances  of  each  gasket  were  experimentally  determined 
(reference  1,  p.  18}.  Each  ring  had  a signal  pickup  feed- 
through which  protruded  through  the  chamber  endplate.  These 
signals  were  fed  directly  down  the  long  50ft  cables  to  the 
instrument  vans  where  they  passed  through  cable  equalizers, 
delays,  attenuators  and  terminations  to  the  oscilloscopes. 

All  intervening  components  were  of  course  properly  matched 
so  that,  from  the  point  of  view  of  the  sensor  input,  the 
effective  load  was  simply  50ft. 
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The  flush-mounted  probes  for  measuring  the  radial  current 
density  are  discussed  in  Appendix  A.  Figure  5-B  is  a photograph 
of  one  of  the  sensors.  The  output  was  tapped  by  direct  connection 
of  a 5 on  cable  to  two  screw  terminals  on  the  back  of  the  sensor 
mount  itself.  This  signal  was  fed  down  the  cable  to  the 
instrument  van  where  it  passed  through  components  similar  to 
those  discussed  in  the  previous  paragraph  for  the  current 
collection  sensors.  The  effective  load  on  each  flush-mounted 
probe  was  50ft. 
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Kamsn  Sciences  Corporation 
ATTN:  Donald  H.  Bryce 
ATTN:  John  R.  Hoffman 
ATTN:  Library 
ATTN:  W.  Foster  Rich 
ATTN:  R.  J.  Fisher 
ATTN:  R.  L.  Hartman 
ATTN:  Walter  E.  Ware 
ATTN:  Albert  P.  Bridges 

Litton  Systems,  Inc. 

Guidance  & Control  Systems  Division 
ATTN:  R.  W.  Maughmer 
ATTN:  John  P.  Retzler 
ATTN:  ValJ.  Ashby,  M.S.  67 

Lockheed  Misuileg  ft  Space  Company 

ATTN:  Clarence  F.  Kool,  Dept.  52-11 
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Lockheed  Missiles  6 Space  Co.t  Inc. 

ATTN:  Samuel  I.  Taimuty,  Dept.  85-85 
ATTN:  George  E.  Heath,  Dept.  81-14 
ATTN:  Philip  J.  Hart,  Dept.  81-14 
ATTN:  Edwin  A.  Smith,  Dept.  85-85 
ATTN:  Benjamin  T.  Kimura,  Dept.  81-14 

LTV  Aerospace  Corporation 
Vought  Systems  Division 

ATTN-  Charles  H.  Coleman 

LTV  Aerospace  Corporation 
Michigan  Division 

ATTN:  Tech.  Lib. 

ATTN:  James  F.  8anson,  B-2 

M.I.T.  Lincoln  Laboratory 
ATTN:  Jean  L.  Ryan 
ATTN:  Leona  Loughiin,  Librarian,  A-082 

Martin  Marietta  Aerospace 
Orlando  Division 

ATTN:  William  W.  Mras,  MP-413 
ATTN:  MonaC.  Griffith,  i ib. , MP-30 
ATTN:  Jack  M.  Ashford,  MP-537 

Martin  Marietta  Corporation 
Denver  Division 

ATTN:  Ben  T.  Graham,  M.S.  PO-454 
ATTN:  J.  E.  Goodwin,  Mail  0452 

Maxwell  Laboratories,  Inc. 

ATTN:  Victor  Fargo 

McDonnell  Douglas  Corporation 
ATTN:  Chester  G.  Polak 

McDonnell  Douglas  Corporation 
ATTN:  Stanley  Schneider 
ATTN:  Paul  H.  Duncan,  Jr. 

Mission  Research  Corporation 
ATTN:  Daniel  F.  Higgins 
ATTN:  Roger  Slettner 
ATTN:  Conrad  L.  Longmire 
ATTN:  William  C.  Hart 

Mission  Research  Corporation 
ATTN:  David  E.  Merewether 
ATTN.  Larry  D.  Scott 

Mission  Research  Corporation,  8an  Diego 
ATTN:  V.  A.  J.  Van  Lint 

The  Mitre  Corporation 

ATTN:  Theodore  Jarvis 

Motorola,  Inc. 

Government  Electronics  Division 

ATTN:  James  R.  Black,  M.S.  A-112 
ATTN:  A.  J.  Kordalewski,  Tech.  Info.  Ctr. 

Northrop  Corporation 
Electronic  Division 

ATTN:  John  M.  Reynolds 
ATTN:  Vincent  R.  DeMartlno 
ATTN:  Boyce  T.  Ahlport 
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Northrop  Corporation 

ATTN:  Orlie  L.  Curtis,  Jr. 

ATTN:  David  N.  Poeock 
ATTN:  Library 

Northrop  Corporation 
Electronic  Division 

ATTN:  Joseph  D.  Russo 

Physics  International  Company 

ATTN:  Doc.  Con.  for  John  H.  Huntington 
ATTN:  Doc.  Con.  for  Philip  W.  Spence 
ATTN:  Doc.  Con.  forlauD.  Smith 
ATTN;  Doc.  Con.  for  Charles  H.  Stallings 

Pulsar  Associates,  Inc. 

ATTN:  Carleton  H.  Jones,  Jr. 

R & D Associates 

ATTN : Leonard  Schlessluger 
ATTN:  S.  Clay  Rogers 
ATTN:  William  R.  Graham,  Jr. 

ATTN:  Richard  R.  Schaefer 
ATTN:  William  J.  Karzas 

The  Rand  Corporation 
ATTN;  Culien  Crain 

Raytheon  Company 

ATTN:  Gajanan  H.  Josbi,  Radar  Sys.  Lab. 

Raytheon  Company 

ATTN:  Harold  L.  Fieachcr 
ATTN:  Jame.  R.  Wcckback 

RCA  Corporation 

Government  & Commercial  Systems 
ATTN:  George  J.  Brucker 

RCA  Corporation 

ATTN:  K.  H.  Zalningcr 

RCA  Corporation 

Government  & Commercial  Systems 
ATTN:  Eleanor  K.  Daly 
ATTN;  Andrew  L.  Warren 

RCA  Corporation 

ATTN:  E.  Van  Keuren,  13-5-2 

Research  Triangle  Institute 

ATTN:  8ec.  Officer  for  Eng.  Dlv.,  Mayr ant  Simons,  Jr. 

Rockwell  International  Corporation 
A'lTN:  Donald  J.  Stevens,  FA-70 
ATTN;  L.  H.  Pinson,  FB-41 
ATTN:  N.  J.  Rudie,  FA-53 
ATTN:  J.  Spetz 

ATTN:  George  C.  Messenger,  FB-61 
ATTN:  James  E.  Bell,  HA-10 

Texas  Instruments,  Inc. 

ATTN:  Donald  J.  Manus,  M.S.  72 

Texas  Tech  University 

ATTN:  Travis  L.  Simpson 


Rockwell  International  Corporation 
ATTN:  John  F.  Roberts 

Rockwell  International  Corporation 
ATTN:  T.  B.  Yates 

Sanders  Associates,  inc. 

ATTN:  James  L.  Burrows 
ATTN:  Moe  L.  Altel,  NCA  1-3236 
ATTN:  K.  G.  Despathy,  Sr.,  PE,  1-6270 

Science  Applications,  Inc. 

ATTN:  William  L.  Chadsev 

Science  Applications,  Inc. 

Huntsville  Division 

ATTN:  Noel  R.  Byrn 

Science  Applications,  Inc. 

ATTN:  j;  Roger  Hill 

Science  Applications,  Inc. 

ATTN;  Charles  Stevens 

Sidney  Frankel  & Associates 
ATTN : Sidney  Frankel 

Simulation  Physics,  Inc. 

ATTN:  Roger  G.  Little 

The  Singer  Company 

ATTN:  Irwin  Goldman,  Eng.  Management 

Sperry  Flight  Systems  Division 
Sperry  Rand  Corporation 

ATTN:  D.  Andrew  Sc  how 

Sperry  Rand  Corporation 
Unlvac  Division 

ATTN:  James  A.  Inds,  M.S.  41T25 

Sperry  Rand  Corporation 
Sperry  Division 

ATTN:  Paul  Marraffino 
ATTN:  Charles  L.  Craig,  EV 

Stanford  Research  Institute 

ATTN:  Robert  A.  A r ml  stead 
ATTN:  Setsuo  Dalrtkl 
ATTN:  Mel  Bernstein 
ATTN:  Philip  J.  Dolan 
ATTN:  Arthur  Lee  Whitson 

Stanford  Research  Institute 

ATTN:  MacPherson  Morgan 

Sundstrand  Corporation 

ATTN:  Curtis  B.  White 

Systems,  Science  & Software 
ATTN:  David  A.  Meskan 

Systems,  Science  t Software,  Inc. 

ATTN:  Ira  Katz 
ATTN:  Andrew  R.  Wilson 

Systron-Donner  Corporation 
ATTN:  Harold  D.  Morris 
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TRW  Systems  Group 

ATTN:  Robert  M.  Webb,  M.S.  Rl-1150 
ATTN:  Tech.  Info.  Ctr.,  S-1930 
ATTN:  William  H.  Robinette,  Jr. 

ATTN:  Jerry  I.  Lubell 
ATTN:  Aaron  H.  Narevsky,  Rl-2144 
ATTN:  Lillian  D.  Singletary,  Rl-1070 
ATTN:  Philip  R.  Garder,  Rl-1028 
ATTN:  Richard  H.  Klngsland,  Rl-2154 
ATTN:  Donald  Jortner 
ATTN;  Paul  Molmud,  SI- 1196 
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TRW  Systems  Group 

ATTN:  Donald  W.  Pugsley 

United  Technologies  Corporation 
Norden  Division 

ATTN:  Conrad  Corda 

Westlnghouse  Electric  Corporation 

ATTN:  Henry  P.  Kalapaca,  M.S.  3526 


TRW  Systems  Gro-^p 
Son  Bernardino  Operations 
ATTN:  J.  M.  Gorman 
ATTN:  Earl  W.  Allen 
ATTN:  JohnE.  Dahnke 


